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ABSTRACT 

We examine the fate of fast electrons (with energies E > 10 eV) in a thermal gas of primor- 
dial composition. To follow their interactions with the background gas, we construct a Monte 
Carlo model that includes: (1) electron-electron scattering (which transforms the electron ki- 
netic energy into heat), (2) collisional ionization of hydrogen and helium (which produces 
secondary electrons that themselves scatter through the medium), and (3) collisional exci- 
tation (which produces secondary photons, whose fates we also follow approximately). For 
the last process, we explicitly include all transitions to upper levels n < 4, together with 
a well-motivated extrapolation to higher levels. In all cases, we use recent calculated cross- 
sections at E < 1 keV and the Bethe approximation to extrapolate to higher energies. We 
compute the fractions of energy deposited as heat, ionization (tracking HI and the helium 
species separately), and excitation (tracking HI Lya separately) under a broad range of con- 
ditions appropriate to the intergalactic medium. The energy deposition fractions depend on 
both the background ionized fraction and the electron energy but are nearly independent of 
the background density. We find good agreement with some, but not all, previous calculations 
at high energies. Electronic tables of our results are available on request. 
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1 INTRODUCTION 



Fast electrons, typically generated by high-energy photons or cos- 
I ■ mic ray collisions, are crucial for a wide range of astrophysical 
' problems. For example, cosmic ray io nization is impo rtant for the 
thermal balance of interstellar clouds ( Spi tzer & Scotti 1969). elec- 
T* ■ trons produced by the Comptonizatio n of 7-ray photo ns affect the 
. ear ly history of supernova remnants dXu et al.ll 199 lh . and broad- 
, line emission regions around quasars are likely ionized by hard 
5^ ' photons from those sources, ma king the fate of the liberated elec- 
trons of paramount importance ( Shu ll & van Steenberdl 19851) . 

As such, there have been numerous studies of the inter- 
actions between these electrons and a background gas, through 
the processes of collisional excitation, ionization, and electron- 
electron sca ttering. These calculations have included both analytic 
approaches (Spitzer & Tomasko 1968; Spitzer & Scott 1969; Jura 
ll97ll ; lBergeron & Collin-Souffrinlll973l;IXu & McCrav 1991) and 
numerical explorations dHabing & Goldsmith] 1 197 it IShullI 1 19791 ; 
IShull & van Steenberd [l985l ; IValdes & Ferrarall2008l) . Most have 
focused on interactions with atomic or ionic gases, relevant espe- 
cially to low-density material in the interstellar (or intergalactic) 
medi a. Others have consid ered the additional effects of molecules 
(e.g.. [DalgarnoetalJl999h . 

Recently, the fate of X-rays in the high-redshift intergalac- 
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tic medium (IGM) has become an important question. Before the 
reionization of HI, ultraviolet ionizing photons are trapped near 
their sources, but X-rays can travel much larger distances through 
the IGM. As such, they are thought to provide the most important 
radiation background for the bulk of the IGM, slowly ionizing the 
mostly neutral gas far from stars. However, most of their energy is 
deposited indirectly, through the fast electrons generated by photo- 
ionization. For example, the electrons produced by X-rays from the 
first hard so urces (quasars, supernova remnan ts, or even Popula- 
tion III stars; Oh 2001;; IVenkatesan et alJuOOlf) are most likely re- 
sponsible for heating the IGM to T ~ 10 00 K before reionization 
dKuhlen & Madaull2005l ; lFurlanettdl2006t) . 

This heating and i onization could influence future genera- 
tions of structure (e.g.. iRicotti et all |2002| ; bh & Haimaij|2003l) 
and it has important observational consequences. For example, 
it can dramatically c hange the highly-reds h ifted 21 cm signal 
from the early IGM jFurlanetto et alj l2006t iKuhlen et al] l2006t 
IPritchard & Furlanettol 120070 . and it could even provide a sig- 
nature of exotic physics during the Dark Age s before the first 
sourc e s appear dChen & Kamionkowskil |2004|; iFurlanetto et al.1 
120061 ; iMapelli et alj 120061 ; IValdes & Ferrard l2008h . Moreover. 
HI Lya photons produced as the secondary electrons colli- 
sionally excite atoms in the background gas can greatly af- 
fect the 21 c m signal by changing the hyperfine level popula- 
tions of HI dChuzhov et alj [20061 ; IPritchard & Furlanettol 120071 ; 
Chen & Miralda-Escude 2008), and the secondary ionizations may 
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play an important role in t he HI (and Hel) reionization of the IGM 
dRicotti & Ostrikej|2004l ; iRicotti et al]|2005l ; IVolonteri & Gnedinl 
l2009h . 

Here we revisit this problem using a Monte Carlo model (de- 
scribed in Sj2]l that incorporates the most recent cross-sections for 
interactions between X-rays and primordial gas. We review the rel- 
evant physics of collisional ionization, excitation, and electron scat- 
tering in ^3l5l and comment on some of our approximations in ij6] 
We present our principal results, together with a comparison to pre- 
vious work, in S|7]and a simple example of their utility for heating 
of the high-redshift IGM in Sj8] Finally, we conclude in ^9] 



2 THE MONTE CARLO MODEL 

Because of the wide range of interactions available to fast elec- 
trons, we use a Monte C arlo model to track the i r fates . Our pro- 
cedure is v ery similar to Shull & van Steenberd jl985h (see also 
IShullI 19791 and lValdes & Ferraral2008h . except that we use updated 
cross-sections and include (and track) more interaction processes. 
An a lternative, fully analytic, approach uses the d egradation equa- 
tion dSpencer&Fanolll954 IXu & McCravlll99ll) , but we find that 
to be more cumbersome. 

As inputs, the model requires only the number densities and 
ionized fractions of hydrogen and helium. We will see that the ion- 
ized fractions have an important effect on the results but that the 
absolute number densities have only a slight impact (through the 
Coulomb logarithm). For simplicity, we will typically assume that 
the HI and Hel fractions are identical (Xi) and that the Hell fraction 
is 1 — Xi. The presence of Helll does not significantly affect our re- 
sults. As a fiducial value, we use the mean density of the Universe 
at z = 10 (according to the cosmological parameters recommended 
bv lDunklev et alj 2009), but taking any density n < 10 s cm 3 only 
changes our reported values by a few percent. Note that our model 
assumes a static IGM and does not self-consistently account for the 
ionization (and heating) produced by each electron. 

The model begins with a single electron of energy E > 
10.2 eV; lower energy particles can only interact with the elec- 
tron gas and so automatically deposit all of their energy as heat. 
Our goal is to compute the fraction of this energy that is deposited 
in ionizing each of the relevant IGM species (HI, Hel, or Hell), 
the fraction of energy that is lost to collisional excitations produc- 
ing photons with E < 13.6 eV (and additionally the fraction that 
ends up in HI Lya photons), and the fraction of energy that is de- 
posited as heat in the IGM. We also follow collisional excitations 
that produce higher energy photons (from helium), but the result- 
ing photons are re-absorbed by the IGM and so are not the ultimate 
products of the process. 

To this end, we compute the cross-sections (weighted by the 
number densities of each species) for the electron to interact with 
the IGM in several different ways: collisional ionization (of any 
of the three species above), collisional excitation (again of any of 
the three species, and including explicitly all levels n < 4), and 
electron-electron collisions. Below we describe our treatment of 
each of these processes in detail. From these cross-sections we ran- 
domly choose one such process for the electron to undergo and 
update its energy, and the energy deposition fractions, accordingly. 
If it ionizes an atom or ion, we add a new secondary electron of the 
appropriate energy to the array of particles; if it collisionally excites 



helium, we add a new photon to the array0 We repeat this process 
until the primary electron's energy falls below 10.2 eV; beyond that 
point, the only available process is a collision with an electron, so 
we assign all its energy to heat. 

We then track each of the secondary electrons through the 
same machinery, and finally we randomly determine (from the 
photoionization cross-sections weighted by number density) the 
species that each secondary photon ionizes and follow the result- 
ing electrons through their entire energy loss cascade. At each step, 
we track the energy lost to each of the aforementioned processes. 

For our final results, we follow 10 J input electrons at each of 
258 energies (logarithmically spaced from 10 eV to 9900 eV) at 
each of fourteen ionized fractions from Xi = 10~ 4 to 0.999. This 
encompasses the range expected in the early IGM (where x% is set 
by the relic density following recombination; Sea ger et al.l 19991) up 
to the point at which the neutrals are no longer significant energy 
sinks. 

To generate random numbers, we use the Mersenne 
Twister algorithm, which has a period of at least 2 19937 — 1 
dMatsumoto & Nishimurj|l998f) . 



3 COLLISIONAL IONIZATION 

For E < 1 keV, we take the collisional ionization cross-sections 
from the CCC database^ an online collection of cross-sections 
calculated with the convergent close-coupling (CCC) method. The 
CCC approach is accurate whenever the target particle can be well- 
modelled by one or two valence electrons above a Hartree-Fock 
core; obviously this is an excellent approximation for hydrogen and 
helium. The relevant physics, and references to many of the orig- 
inal p apers for hydrogen and helium, can be found in iBrav et al.l 
1 120021) . We use cubic splines to interpolate the database values. 

At higher energies, w e assume th at cross-sections follow the 
Bethe approximation limit (Beth el 19301) . Then the cross-section for 
ionization from the ground state in species i is 



a x ~ ¥k ( At In — + B, 
E V 2Ei 



2 

Tva , 



(1) 



where ao is the Bohr radius and A4 and Bi are coefficients. These 
are usually fixed by demanding that the cross-section map onto the 
first Born approximation at high energies. This asymptotic behavior 
compares well to the calculated behavior at E ~ 1 keV, and we ex- 
trapolate to higher energies by fitting a function of this form to the 
uppermost energy bin in the CCC results. The resu lting parameters 
Ah and Bh differ from the analytic estimates of Jo hnsonl d 19721) 
by ~ 15% and 30%, respectively, but their exact values make lit- 
tle difference to our final results (largely because it is the relative 
importance of each ionization and excitation process that matters). 

Whenever a species is collisionally ionized, it also produces a 
secondary electron, so we must select the final energies of the inci- 
dent and ejected elect rons. We use the probability distributions of 
iDalgamo etafl h999T) . which are adapted from the measurements 
of lOpaletalJ<197ll) : see also lShullI J 19791) for a discussion of the 
energy spectra of secondaries. This prescription assumes that the 
probability for the secondary to have an energy e is proportional to 



Pie) 



l+te/ei) 2 - 1 ' 



(2) 



1 Collisional excitation of hydrogen produces a photon with E < 13.6 eV 
that no longer affects the IGM. 

2 See http://at 6m.curtin.edu.au/CCC- WWW/ 
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where e, = 8, 15.8, and 32.6 eV for HI, Hel, and Hell, respec- 
tively. (We always identify the "secondary" as having the lower 
energy of the two final photons, so that e < (E — Ei)/2.) Note 
that secondary electrons are typically ejected with a modest en- 
ergy somewhat below the ionization threshold of their original host 
atom: the median energies are 7.2, 14.2, and 28.5 eV for the three 
species, alth ough there is a tail to much higher e. As discussed in 
IShullldl979h . the mean secondary energy increases logarithmically 
with the incid e nt ene rgy for fast collisions^ 

Ijohnsonl d 1972b provides an alternate estimate of the sec- 
ondary electron energy distribution for HI by constructing analytic 
approximations to the HI ionization and excitation cross-sections. 
In his calculation, the leading order behavior is oc (1 + E / Ei)~ 2 , 
although he also includes several other correctiv e terms propor - 
tional to higher powers of this same factor (see also lOmidvad 1969b . 
This has a similar (though not identical) shape to our choice, and 
there is no simple way to extend it to the other relevant ions. 



4 COLLISIONAL EXCITATION 

We include collisional excitation processes in a similar way to col- 
lisional ionization, by making use of the CCC database. This pro- 
vides cross-sections for all excitations from the n = 1 to n = 2, 3, 
and 4 states of HI, Hel, and Hell, including separately all angular 
momentum sublevels as well as separate estimates for the singlet 
and triplet configurations in the case of Hel. We again extrapolate 
to E > 1 keV using the Bethe approximation. 

This explicit separation is useful to us because some ap- 
plications (particularly estimates of the high-redshift 21 cm sig- 
nal) require the production rate of HI Lya photons through col- 
lisio nal excitation. For the n < 4 l evels, this i s relatively easy 
(see IPritchard & Furlanettoll2006l and|Hirata 2006: for discussions 
of analogous photo-excitation processes)Q Atoms excited to the 2s 
state obviously produce no Lya photons. Atoms excited to the 3s or 
3d levels must either decay through the 2p state (producing a Lya 
photon) or directly to the ground state. But in the latter case the 
resulting line photon will quickly be re-absorbed by a nearby atom 
(assuming that neutrals are encountered before the photon redshifts 
out of resonance), and it will eventually cascade into a Lya photon. 
On the other hand, atoms excited to the 3p state must decay to 2s 
and then to the ground state via two-photon decay. Thus they do 
not produce Lya photons at all. 

The cascade possibilities become more complicated at n = 4. 
For example, excitation to 4s can result in decay directly to the 
ground state (which is irrelevant for our purposes, since the pho- 
ton is re-absorbed) or indirectly through the 3p state (in which case 
no Lya photon is produced) or the 2p state (which produces a Lya 
photon). The relative fractions of these are given by the decay prob- 
abilities, 



PLya(4s) = 



-4 



4s — 2p 



A'l 



(3) 



_l_ Ae ' 
4s — 2p ' A 4s-3p 

where A e rlm is the spontaneous emission coefficient from level n to 
m. Similar exercises yield the decay probabilities for the other n = 
4 sublevels. It is important to include excitations to the different 
n > 2 angular momentum sublevels separately, because they feed 



3 Note that the exponent of 2.1 and coefficients ei in equation (2) differ 
slightly from the fits of lShuJ ll979l) . 

4 Note that the discussion below ignores the possibility of collisional de- 
excitation and so will break down in sufficiently dense gas. 



into Lya differe ntly; most previou s work has included only n < 2 
(e.g., Shull & van Steenberg 1985), averaged over all the sublevels 
(e.g.. lXu & McCravl 1 199 lb. or include d only excitation to the np 
sublevels (e.gJValdes & Fer rar3l2008b . 

Note that, in the case of Hel and Hell excitations, the decay 
photons can ionize other elements. In detail, we should follow the 
radiative cascades of these excited atoms or ions, propagate each 
decay photon through the IGM, and follow the resulting ioniza- 
tions. However, we take a cruder approach and simply assume that 
each excited atom or ion immediately decays to the ground state, 
producing a photon equal to the energy of the excitation. We then 
assume that this photon ionizes another atom and follow the result- 
ing photo-electron. We use the photoionization cross-sections from 
IVerner et al.l j 1996b for this purpose. 

We therefore eventually recycle all the Hel and Hell excitation 
energy into ionization or heating, rather than allowing a fraction to 
escape. However, we find that < 1.5% of the initial electron energy 
is typically lost to Hel excitation (and even smaller than that for 
Hell) in the high energy limit, and most of that is direct excitation 
to 2p, for which our assumption is accurate. This approximation 
is therefore a small correction to our results (comparable to many 
other processes that we ignore, and smaller than uncertainties in 
our cross-sections). It will, however, affect the results more strongl y 
near the Hel threshold (e.g., Fig. 1 o f lShull & van Steenberg! 1985). 



4.1 Excitation to n > 4 states 

The CCC database does not include cross-sections for transitions 
to n > 4 states, and in any case it is impractical to include arbi- 
trarily many final states. We therefore approxima te the effects of 
highe r levels based on the analytic calculations o f IJohnsonl fl972h 
and lXu & McCravl dl99lh - 

In the Bethe approximation, the leading coefficient A n (for 
excitation from the ground state to the nth level) is proportional to 
the oscillator strength of the equivalent radiative absorption transi- 
tion, /„. This depends on the quantum number of the upper level 
through /„ cx g n /{'nk n f', where k n — 1 — (1/n) 2 and g n is a 
correction factor (the "Gaunt factor") that depends on the quantum 
mechanical details of the interaction. Thus we can approximate the 
total cross-section to all levels above n — 4, relative to that for 
n = 4, as 



^=60V(fc„n)- 3 

n — 5 



1.58. 



(4) 



Because of the details that we have ignored (namely the 
Gaunt factor), the actual cross- section of the higher l evels is some- 
what smaller than this estimate. IXu& McCravl jl 99 lh explicitly in- 
cluded all n < 10 in t heir calculation (utilizing the analytic approx- 
imate cross-section of Ijohnsonll 1972b : they found that the n — 5- 
10 states increased the interaction rate by a factor of 1.17, as op- 
posed to the 1.29 suggested by truncating the sum in equation {4]l 
at n = 10. We therefore increase the n — 4 cross-sections by a 
factor (1.17/1.29) x 1.58 = 1.43 for a fiducial estimate of the 
importance of higher-level transitions. Fortunately, varying this en- 
hancement level by a factor of two affects our heating, ionization, 
and excitation fractions by less than a percent. 

We further assume that the fraction of such excitations produc- 
ing Lya p hotons is identical to that for n — 4, which is a reasonable 
estimate dPritchard & Furlanettdl2006l) but does not substantially 
affect our results anyway. 

We also apply this same n > 4 enhancement to Hel and Hell, 
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although this is no more than a guess. Again, this makes only a 
negligible difference to our results. 



5 ELECTRON-ELECTRON COLLISIONS AND 
HEATING 

A fast electron will share its kinetic energy with its surround- 
in gs by scattering off o f ambient electrons. The energy loss rate 
is ispitzer & Scottll 19691) 

dE _ 47re 4 n e lnA 

7, 7 

at m e w 



where In A is the Coulomb logarithm (to be discussed in §63]be- 
low) and w is the velocity of the electron. However, this near- 
continuous slowing-down (thanks to the long-range interactions of 
the Coulomb force) must be discretized to be included in our Monte 

Carlo formalism. 

We follow lHabing & Goldsmith! Jl97 lh and IShullI d 19791) by 
casting this process in the following form. First, we assume that 
the electron loses a fixed fraction / of its energy in each interaction. 
We then rewrite the energy loss rate as 

dE ^ fE 
dt t co n 

where the collision time is t co \\ — 1/ (wn e a ee ) and we have let 
the interaction cross-section be tj ee . Rearranging, and taking the 
non-relativistic limit E — mw 2 /2, we obtain 



fE x wn e a e 



(6) 



407re 4 



In A 



0.05 
/ 



(7) 



Although the parameter / is artificial, it does not affect our 
calculations: all models with / < 0.05 yield identical results, 
within the errors expected in our Monte Carlo simulations. 



6 APPROXIMATIONS 

Before proceeding to our results, we now discuss a few sub- 
tleties that arise in the calculation. In addition to the processes 
that we discuss in detail below, our model also ignores several 
subdominant mechanisms that have been studied before; we list 
them here for completeness. First, we ignore the double ioniza- 
tion of neutral helium, which occu rs at a rate ~ 2% that of sin- 
gle ionization jDalgarnoetZTl 999). We also ignore Coulomb col- 
lisions with protons and recombinations. IValdes & Fe rrara (2008) 
included these two latter processes and showed that they are negli- 
gible. 



ti oss ,H ~ 5 x lO 5 ^ 1 ( — 



where we have used the Coulomb logarithm appropriate to the IGM 
at z = 10. Clearly, even with an ionized fraction Xi > 10~ 3 , this 
timescale is only a small fraction of the Hubble time, so the primary 
photoelectron will rapidly lose its energy once it is created. 

At smaller ionized fractions, electron scattering may be slow. 
However, with the Bethe approximation form for the collisional 
ionization cross-section, and assuming that AE ~ Eh in each 
interaction, the energy loss timescale from ionization is 

E \ 3/2 (\ + z\- 3 

( — ) yf ' (9) 
where xh = 1 — Xi and we have evaluated the logarithmic factor at 
E = 1 keV as well. The two processes are comparable when Xi ~ 
0.01. Thus ti oss < 5 x 10 5 yr for keV photons: much smaller than 
the Hubble time, so our instantaneous approximation is reasonable. 
The only exceptions are very high energy electrons (which result 
only from high energy X-rays, to which the Universe is transparent 
anyway). 

Of course, if the underlying ionization fraction changes 
rapidly even this timescale may not be short enough. While such 
changes are unlikely to occur faster than the Hubble time on a 
global scale, the ionized fraction can certainly change rapidly on 
a local level when a bright source ionizes its environs. We there- 
fore urge caution in using our results under such circumstances. 

6.2 Inverse Compton scattering 

We have also ignored energy losses due to inverse Compton scat- 
tering from CMB photons. The emitted power due to this process 
is dE/dt — — (4/3)oTC?7rad7 2 /3 2 , where <tt is the Thomson 
cross-section, (7 ra d = glTcmb i s the energy density of the CMB, 
P = w/c, and 7 2 = 1/(1 - g 2 ) (Rvb icki & Lightmaril979h . For a 
non-relativistic electron, the timescale for inverse Compton cooling 
is therefore 



tc 



3m e c 



8cftciTA 



10 ! 



CMB 



yr. 



(10) 



Note that this is independent of energy, so inverse Compton cooling 
is most important for high-energy electrons, where the other loss 
timescales are large. 

Because t comp is typically much longer than the loss timescale 
due to collisional processes (which is bounded from above by eqs.[8] 
and [9}, and because of the strong redshift dependence, we have 
ignored inverse Compton scattering in our calculations. It can be 
roughly included by assuming that a fraction ti oss /t corap of the 
electron energy is lost to the CMB and reducing the other frac- 
tions accordingly (note that their proportions will not change sig- 
nificantly, unless Xi <C 10~ 2 and the energy is near a line thresh- 
old). 



6.1 Timescales 



6.3 The IGM density and temperature 



Our model assumes that all of the photon's energy is deposited in 
the IGM instantaneously. In situations where that energy deposi- 
tion occurs over a long time period, this may not be a good approx- 
imation (for example, if the density field or ionized fraction evolve 
rapidly compared to the energy deposition timescale). We can es- 
timate the relevant timescales by considering electron-electron in- 
teractions and collisional ionization of HI, which are typically the 
most important processes. For the former, we find that 



E 
dE/dt 



5 x lO^av 



(—) 



3/2 



yr, (8) 



In all of our calculations, we assume that the background gas has 
a density equal to the cosmic mean at z = 10. But, to a good 
approximation, our results are independent of the absolute density, 
because the interaction rates all scale linearly with it. However, the 
Coulomb logarithm in equation Q does have an implicit density 
dependence that breaks this convenient feature. If one considers 
only the scattering of the photoele ctron from dis crete particles, the 
appropriate Coulomb logarithm is lSpit zeJl962l) 



In A 



(fcT) a 

7T7i.ee 6 7e ' 



(ID 
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where n e is the electron density, T is the gas temperature, and In -y e 
is Euler's constant. In this regime, this factor is independent of the 
electron's energy, so a ee oc E~ 2 and heating rapidly becomes less 
important at high energies compared to ionization and excitation. 

However, the photoelectron can also have collective interac- 
tions with the plasma as a whole. This process becomes more im- 
portant at high energies. The effective Coulomb logarithm, includ- 
ing both individual co llisions and these collective interactions, is 
JSchunk&Havs1ll97lh 



In A = In 



2E / 47rn e e 2 X ~ 1/2 



(12) 



so long as E > me i /2h 2 « 13.7 eV0 Here, the factor Ce = 
7.40 x 10 _11 (n e /cm" 3 ) eV. In this case, a ee oc E~ 2 lnE, so 
heating is somewhat more important at high energies, although o ee 
still falls off by one power of energy faster than the ionization and 
excitation cross-sections. Note that, when plasma effects are in- 
cluded, the cross-section is independent of the temperature of the 
gas, so we need not specify it in our calculations. 

We use equation d 1 2b in our fiducial calculations, following 
IXu & McCravl dl99 lh . Ignoring collective effects reduces the frac- 
tion of incident energy deposited as heat at E > 100 eV by a few 
percent, and it correspondingly increases the fractions of energy 
deposited in ionization and excitation. 

The density can also play a second role in collisional de- 
excitation of the target atoms and in populating the upper levels of 
the atoms (provided also that the temperature is sufficiently large). 
This requires n 3> 10 s cm" 3 so is not important for the IGM. 
But we caution readers against trusting our results in such environ- 
ments. 



7 RESULTS 

We now present our main results. As a reminder, we use 10 
Monte Carlo trials for each electron energy (which ranges from 
10-9900 eV), and we explicitly follow all secondary electrons and 
photons. We take / = 0.05 and set the density of the background 
gas to be the mean cosmic density at z — 10, but we saw above 
that the results are very insensitive to these parameters. The pri- 
mary input parameters are then the initial electron energy and the 
ionization fractions. For the latter, we let Xi be the density of HII 
relative to the total hydrogen density and assume that this is also 
the fraction of helium in the form of Hell, with zero Helll. 

We present our results in terms of energy deposition fractions, 
with /i on , /heat, and /excite the fractions of the initial electron en- 
ergy that goes into ionization, heating, and HI line photons gener- 
ated by collisional excitation[f| (Note that line photons from helium 
are transformed into ionizing photons and so deposit their energy in 
other ways, implicit in our code.) We also let /L ya be the fraction 
of energy deposited in HI Lya photons. We have verified explicitly 
that our code conserves energy throughout the entire interaction cy- 
cle. 

Figure[T]shows these fractions, as a function of photon energy, 
in a medium with Xi = 0.01. The thick solid, dashed, and dotted 



5 We are interested in slightly lower energies as well, but the dynamics 
there are almost completely dominated by the structure of the line spectrum 
of HI rather than the details of <r ee . 

6 We caution the reader that our parameters do not include the effects of 
the initial ionization event that generates the electron. 
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Figure 1. Energy deposition fractions as a function of initial electron en- 
ergy, for an IGM with Xi = 0.01. The thick solid curve shows the fraction 
of energy used in ionization; the thin solid and dot-dashed curves show the 
fractions used for ionizing HI and Hel, respectively. The thick dashed curve 
shows the fraction deposited as heat through electron-electron interactions. 
The thick dotted curve shows the total energy lost to HI line photons; the 
thin dotted curve shows the fraction in HI Lya photons. 



curves show /i on , /heat, and /excite, respectively. The thin solid 
and dot-dashed curves show explicitly the fractions of energy going 
into HI and Hel ionization, respectively. (The fraction going into 
Hell ionization is <C 1% throughout.) Finally, the thin dotted curve 
shows f Lya . 

The structure of the curves is relatively easy to understand. 
Electrons with E < 10.2 eV are unable to interact with any atoms 
or ions, so all of the energy is deposited as heat. As E increases, 
more and more excitation and ionization processes become avail- 
able, so the energy injected as heat decreases. Interestingly, even 
with this large of a neutral fraction, the individual line thresholds 
do not introduce discrete features into the deposition fractions, and 
all of the parameters depend smoothly on energy. 

At higher energies, where no additional processes become 
available, the fractions vary only slowly, eventually approaching 
reasonably constant values at E ~ 1-10 keV, where ioniza- 
tio n, heating, and excitation split roughly equally. As pointed out 
by IShull & van Steenberd dl985h . the naive expectation from the 
Bethe approximation (with <Ji, e oc E~ In E for ionization and ex- 
citation and (j ee oc E~ 2 In E for heating), that less and less heating 
should occur at high photon energies, is false. This is because each 
ionization produces a moderate energy secondary electron (typi- 
cally ~ 10 eV, but occasionally much larger). A large fraction of 
the primary's energy is lost through these intermediaries, who in 
turn lose most of their energy to heat, making the behavior at high 
energies much less variable than one would expect. 

Figures [2] and [3] show the energy deposition fractions for a 
range of Xi. The first panel shows /i on . For small Xi, the ionization 
thresholds imprint features on the curves. At high energies, /i on ~ 
0.4 at small ionized fractions and rapidly approaches zero as Xi 
increases. The behavior at very high energies sets an upper limit to 
the total energy invested in ionization. 
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Figure 2. Energy deposition fractions in ionization (left) and heating (right) as a function of initial electron energy. The curves show results for several different 
IGM ionized fractions, as labelled in the left panel. 
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Figure 3. Energy deposition fractions in collisional excitation (left) as a function of initial electron energy, and the fraction of this energy that goes into HI 
Lya photons (right). The curves show results for several different IGM ionized fractions, as labelled in the left panel. Note that we do not show curves for 
Xi = 0.5 or Xi = 0.9 in the right panel, as they are very noisy at low energies (although the high-energy limit is well-behaved). 



The left panel of Figure [5] shows the corresponding fractions 
lost to collisionally excited HI line photons. Note the line struc- 
ture apparent at xi < 10~ 3 ; the features are separated by ~ 10 eV 
and correspond to multiple excitations of HI atoms. Although the 
structure around these features appears noisy, it is real and remains 
unchanged in other Monte Carlo tests. However, as noted above 
these features blend together by Xi ~ 10 -2 . Note that collisional 
excitation becomes important at lower photon energies than does 
ionization, because the lines require less energy; in nearly neutral 



gas, these lines are clearly the dominant process even near the ion- 
ization threshold. 

Interestingly, the right panel of Figure [5] shows that the ratio 
fhya //excite is nearly independent of both energy and Xi: it de- 
pends almost entirely on the atomic physics of HI. The weak energy 
dependence near threshold occurs because the different levels feed 
into Lya in different ways. The (very) weak dependence on Xi, 
with the Lya fraction increasing slightly with Xi, is probably be- 
cause efficient electron-electron interactions in highly-ionized gas 
bring electrons to the Lya threshold energy more rapidly. The over- 
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all fractio n of ~ 80% can easily be es timated from the oscillator 
strengths jPritchard & Furlanettcll2007l) . 

The remaining energy, shown in the right panel of Figure [2] 
heats the IGM. The curves vary only slowly at large E, but the 
behavior at E < 100 eV is complicated by the ionization and line 
thresholds. Note that the features here are caused by the combina- 
tion of these processes. 



7.1 Comparison to previous work 

Calculations similar to ours have appear ed several times over 
the pa s t four decades. The m ost r ecent are Shull & van Steenberd 
dl985l),IXu & McCravl fl99lh . and lValdes & Ferraral j2008tT 

iXu & McCravl dl99lh used the degradation equation 
53, 



( Spencer &Fanoll 19541) . a fully analytic technique, to compute the 
energy deposition fractions. They took the approximate ( but fully 
analyt ic) ionization and excitation cross-sections from IjohnsonI 
dl972l) . and explicitly included all excitations with n < 10, 
although only for a pure hydrogen gas. 

We have generated a similar set of scenarios to lXu & McCrav 
dl99ll) for comparison purposes, in which we include only hydro- 
gen target atoms, use their Coulomb logarithm, and extrapolate the 
excitation cross-sections to only n < 10 levels. Over the range 
Xi = 10 _4 -0.9, we find generally go od agreement betwe en the 
two sets of results at E — 2 keV (where Xu & McCravl 1 99 ll report 
detailed deposition fractions). The deviations in the energy deposi- 
tion fractions are always smaller than a few percent (in absolute 
terms), with our /heat typically slightly larger and our other frac- 
tions slightly smaller; the diffe rences a re larg est at large ionized 
fractions. We note that Dalg arnoetal]dl999l) also found slightly 
larger heating fractions than lXu & McCravl dl99lh . At Xi < 0.01, 
all of our results agree to within 1% (again in absolute terms). 
Given that we use completely different methods and cross-sections, 
we regard this agreement as excellent. There is also good qualita- 
tive agreement in the shapes of our energy deposition curves, and 
especially in the features from excitation and ionization that appear 
at small x j. 

In contrast, Shull & van Steenberd < fl985T) used a Monte Carlo 
method very similar to ours, although we have added several inter- 
actions (namely excitation to higher levels) and updated the cross- 
sections for most of the others (including the Coulomb logarithm 
in the heating component). Despite these differences, our results 
agree to similar (few percent) accuracy in the high-energy limit 
(or specifically at 3 keV, where IShull & van Steenberdl 19851 pro- 
vide detailed results). Figure [4] provides a detailed comparison as 
a function of Xj\ the solid and dashed cu rves show our results and 
those from lShull & van Steenberal dl985l) . respectively^ 

As one mig ht expect given their la c k of d ata for excitations 
to high n states, IShull & van Steenberd dl985l) somewhat under- 
estimate the importance of collisional excitation, which can cause 
a discrepancy of up to ~ 6% with respect to our results at small 
neutral fractions. They overestimate /i on by a similar amount com- 
pared to our calculations, but /heat agrees rather well (although our 
results clearly have more curvature than theirs). Overall, we regard 
this a greement as quite good. 

IValdes & Ferrari ( 2008) also used a Monte Carlo model to 
examine the energy deposition fractions for high-energy electrons, 



7 For /Lya we have used their expression for the fraction of energy de- 
posited in collisional excitation of HI, multiplied by a factor 0.8 as in Fig- 
urefj] 
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Figure 4. Comparison of our results (points) to previ ous fitting formu- 
lae fro m lsiiull & va n Steenberd <1985t> (thick curves) and Valdes & Ferraral 
120081) (thin curves). The dotted curves and open diamonds show /heat- 
The solid and dashed curves (crosses and open squares) show /; on and 
fhv a, respectively. Most are evaluated at E = 3 keV, except for those 
from lValdes & Ferrard |200S1) . who use E = 2 keV. 



with E > 3 keV. They included most of the same processes 
we have, albeit with different cross-sections in most cases. The 
most significant difference is their treatment of collisional exci- 
tation; they took values from IStone et al. ( 2002), who computed 
their cross-sections wit h the scaled plane-wave Born approxima- 
tion. IStone et all d2002l) show explicit comparisons with the CCC 
databa se; the CCC values differ by < 20%. However, most signif- 
icantly IStone et al I d2002h only compute excitation cross-sections 
to the np angular momentum sublevels. Although these transitions 
are typ ically dominant , the ot hers are certainly non-negligible. 

Valde s & Ferraral d2008h present detailed results for the high- 
energy limit and offer alte rnate fitting functions to those from 
Shull & van Steenberal dl985l) ; Figure [4] also includes a compari- 
son to their results (do tted curves). The agreement w ith ours is not 
nearly as good as for lshull & van Steenberd dl985l) . They find a 
much lower rate of heating at moderate and large ionized fractions 
(by up to ~ 20% at Xi ~ 0.01) and correspondingly larger rates 
of excitation and (especially) ionization. The agreement is partic- 
ularly poor at Xi < 1 (although the practical importance of this 
regime is fairly small). The source of this discrepancy is unclear; 
our excitation cros s-sections match closely in the high-energy limit 
dStone et alj|2002h . and our ionization cross-sections are in good 
agreement over the entire ener gy range (e.g., comparing the CCC 
database to lKim & Rudll 19941) . 

Another difference between our results and[Valdes & F erraral 

d2008l) is in the fraction of excitation energy deposited as Lya pho- 
tons. As shown in Figure[3] we find /l vo //excite ~ 0-8 throu ghout 
the entire energy range. In contrast, Valde s & Ferraral d2008l) find 
/Lya/ /excite ~ 0.7. This discrepancy probably results from their 
inclusion only of excitations to the np sublevels; the mixing frac- 
tions for the (rarer) excitations to the other sublevels can be larger, 
especially at n = 3 (where atoms in the 3p state cannot decay 
through Lya, but those in the 3s or 3d states must). 
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Figure 5. Energy deposition fractions in HI ionization (left) and heating (right) as a function of electron energy. The curves in each panel show results for 
several different I GM ionized fraction s, with log Xi = —4, —3, —2, and —1, from top to bottom. The solid curves show our calculations; the dotted curves 
show the fits from lRicotti et alj 1 2002); see eqs. 413I14> in the text. 



Over lo wer energies, the agreem ent with previous results 
is less clear. PValdes & Ferrara ( 2008) e xamin ed only the high- 
energy limit, but Shull & van Steenberd Jl985h studied the entire 
range. They prese nted few numerical results at lower energies, but 
iRicotti et alj J2002f) provided t he following fits to their results (see 
also lVolonteri & Gnedinll 20091) : 



rR 

lion, HI 



rR 

J heat 



— m 

+0.39(1 -x° A1 ) 
+[l-(1.0-a:?- ; 



E > 28 eV 



4 (1 



E > 11 eV, 



(13) 



(14) 



and zero otherwise. The terms without energy dependence describe 
the high-energy limiting values discussed above; the final term in 
each expression was added to account for the energy dependence, 
at least roughly. 

Figure [5] compares these expressions to our detailed results 
over the entire relevant energy range. The left panel shows /km, hi, 
while the right panel shows /heat • The sets of curves take log Xi = 
—4, —3, —2, and —1, from top to bottom, with the solid curves 
showing our results and the dotted curves using the fits. Although 
(for the most part) the fits are reasonably accurate at E > 1 keV0 
they provide a relatively poor match at lower energies. We therefore 
caution against these fits for high-accuracy work. 

The advantage of equations J 1 31 14b is that the energy and Xi 
dependence are separable, which simplifies their application to nu- 
merical work. Unfortunately, Figures [2] and [3] show that the func- 
tional dependence over energy can change quite substantially with 



Xi (even leaving aside the features from the ionization and excita- 
tion thresholds). We suspect this is why a single, separable fitting 
function fails, and we have not attempted to find another form. In- 
stead, we recommend interpolating the exact results. 



8 AN EXAMPLE APPLICATION: DISCRETE SOURCES 
IN THE EARLY UNIVERSE 

As a simple example of the utility of our results, we briefly de- 
scribe here an application to X-ray heating by discrete sources in 
the high-z universe. Suppose a source emitting ultraviolet and X- 
ray photons with a power law spectrum, L v cx v~ a , is embedded in 
the IGM. We will take a = 1.5 for concreteness. The photons will 
gradually be absorbed as they stream through the neutral gas sur- 
rounding the source: ultraviolet photons just above the ionization 
edge have a very short mean free path, but higher-energy photons 
can reach much larger distances. The comoving mean free path of 
an X-ray photon with energy £J 7 is 

3 



A, 



4.9x 



-1/3 



1 + 2 



E~, 



Mpc. 



(15) 



15 J V300eV/ 

Thus, for a point at a finite distance from our source, all pho- 
tons with energies _E 7 < -Emm will be strongly attenuated, while 
those above this threshold will be more or less unaffected. We ap- 
proximate this situation by assuming zero transmission below -B m i n 
and complete transmission above that threshold. We then compute 
the net fraction of absorbed energy that is deposited as heatEl 



/he 



Si JE 



1 n l a i (E J )(E J - Ei)f he 



Y^ils dE^ E 1 a n i <j i (E 1 



(16) 



8 The largest deviation in the high energy limit, of ~ 6% in absolute terms, 
occurs in /h oa t when Xi = 0.0 1. Note, however, that our re sult actually 
agrees with the exact result from lShull & v an Steenbera 1 1985) well (0.32 
for both); the error comes from the fitting function itself. 



9 By integrating over a here, we are implicitly assuming the optically thin 
limit for the surviving photons. If a system is optically thick, it will weight 
high-energy photons more strongly. 
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Figure 6. Energy deposition fractions for the secondary electrons integrated over sample spectra (with a = 1.5), assuming optically thin absorption. The left 
panel shows /i onj Hi, the upper right panel shows /L ya , ar| d the lower right panel shows /h ca t- The s olid, long-dashed, shor t-dashed, and dot-dashed curves 
take logXi = —4, —3, —2, and —1, respectively. The solid squares show the high-energy fit from Shull & van Steenberg ( 1985), with the same ordering 
from top to bottom. 



where the su m is over all specie s i and cr; is the photoionization 
cross-section dVemeret alii 19961) . Here the numerator is the inte- 
gral (over all photon energies £ 7 ) of the ionization rate multiplied 
by the fraction of the secondary electron's energy [E 1 — Et) de- 
posited as heat. The denominator normalizes / to the total energy 
deposition rate. We can compute /i on .Hi and fh ya in a similar fash- 
ion, replacing /heat with /i on ,Hi or fhya- (Note that / still does not 
include the energy deposited by the initial ionization. This is why 
the sum of these fractions is much smaller than unity at low ener- 
gies.) 

The curves in Figure [6] show these two quantities at several 
different ionized fractions (loga^ = —4, —3, —2, and —1, from 
top to bottom in the left hand panel). The solid squares show the 
asymp totic high-energy estimates from the lShull & van Steenberd 
( 1985) fitting formulae. (These are the terms in eqs. ll3ll4l without 
any energy dependence.) 

These fitting formulae provide reasonable order-of-magnitude 
estimates so long as E m i n > 300 eV (or slightly higher if Xi ~ 
0.1), although they systematically underestimate /L yQ and /i on ,Hi- 
These deviations are most significant at small ionized fractions, 
where heating is least important. We therefore again recommend 
interpolation of the exact results for high-accuracy work, especially 
if the heating by soft X-rays E < 300 eV is included. 



9 DISCUSSION 

Using a Monte Carlo model, we have re-examined the fate of fast 
electrons scattering through a background gas of primordial ori- 
gin. We included electron-electron scattering as well as collisional 
ionization and excitation of HI, Hel, and Hell, explicitly track- 
ing all levels up to n — 4 and using an analytic extrapolation 
to higher levels. We separately followed all excitations produc- 
ing HI Lya photons, which can be important in modeling the ob- 



servable properties of the IGM at high redshifts (see Kuhl en et al.l 
2006 and Furlanetto et al .120061 for example) and have not been ex- 
plicitly tracked previously except at the highest energies. We used 
recent calculations of ionization and excitation cross-sections at 
E < 1 keV and extrapolated to higher energies using the Bethe 
approximation. 

In highly neutral gas (xi < 10 -3 ), we found that ~ 20% of 
the electron energy is deposited as heat, with the remainder split 
roughly equally between ionization and excitation. In this regime, 
the results are not strongly sensitive to Xi, at least at high ener- 
gies, because most of the heating comes from secondary electrons, 
with energies below 10 eV. At higher Xi, the heating fraction rises 
rapidly, exceeding ~ 65% by Xi ~ 0.1. We find that the excitation 
and ionization energy deposition rates are always comparable, and 
that ~ 80% of the excitation energy goes into HI Lya regardless of 
electron energy and x%. 

Although our calculations used parameters appropriate to the 
low-density IGM, our results may also be applied to denser sys- 
tems, because the density only enters through the Coulomb log- 
arithm affecting the electron-electron scattering rate. Varying the 
density of target atoms by many orders of magnitude only affects 
the energy deposition fractions by a few percent in absolute terms. 
We also note that, when collective plasma effects are included, the 
background temperature be comes irrelevant for the electron en- 
ergies under consideration ( Schu nk & Havsl Il97ll : IXu & McCravl 

[mil). 

For the most part, our results agr ee with previou s estim ates 
from lShull & van Steenber j dl985h and lXu & McCravl Jl99lh . al- 
though we have found some discrepancies with the commonly-used 
fitting formulae from the former. In general, at high energies their 
results slightly underestimate the importance of collisional excita- 
tion but are otherwise accurate. At lower energies, the differences 
in cross-sections become more important and the disc repancies in- 
crease (at least with reference to the fitting formulae of Rico tti et al.l 
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I2002I). How ever, our results show substantial differences with those 
of lValdes&Ferraral (2008), who examined the high-energy limit 
with a Monte Carlo model similar to ours. These discrepancies are 
especially large at moderate and high ionized fractions, where we 
find substantially more heating and less ionization and excitation. 
We also find that a higher fraction of excitation energy is deposited 
in HI Lya photons. The latter is probably due to our different treat- 
ments of collisional excitation (in particular, their neglect of exci- 
tations to states other than the np sublevels), but the source of the 
former is unclear. 

In any case, we advocate interpolation of the exact results 
when high accuracy is necessary, especially because the energy 
dependence is quite significant. Our detailed numerical results are 
available upon request, including tables for the energy deposition 
fractions at a variety of ionized fractions and C code to interpolate 
to arbitrary ionized fractions and electron energies. 
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